Memory Effects in Spontaneous Emission Processes 
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We consider a quantum-mechanical analysis of spontaneous emission in terms of an effective two- 
level system with a vacuum decay rate To and transition angular frequency lua- Our analysis is 
in principle exact, even though presented as a numerical solution of the time-evolution including 
memory effects. The results so obtained are confronted with previous discussions in the literature. 
In terms of the dimensionless lifetime r = tTo of spontaneous emission, we obtain deviations from 
exponential decay of the form 0(1/t) for the decay amplitude as well as the previously obtained 
asymptotic behaviors of the form 0(1/ t 2 ) or 0(1/ t In 2 r) for t>1. The actual asymptotic behavior 
depends on the adopted regularization procedure as well as on the physical parameters at hand. We 
show that for any reasonable range of r and for a sufficiently large value of the required angular 
frequency cut-off ui c of the electro-magnetic fluctuations, i.e. uj c 3> uja, one obtains either a 0(1/t) 
or a 0(1/t 2 ) dependence. In the presence of physical boundaries, which can change the decay 
rate with many orders of magnitude, the conclusions remains the same after a suitable rescaling of 
parameters. 

PACS numbers: 03.65.-w, 12.20.Ds, 32.10.-f 



I. INTRODUCTION 



The development concerning the manipulation of sin- 
gle atoms and their interaction with the electro-magnetic 
field has reached an impressive state of art in recent years 
(see e.g. Refs. [H-Q)- Experimental studies have e.g. 
shown that artificial atoms can lead to a Lamb shift of 
the order of a few per cent of the typical emission line 
[H ■ Rather old issues concerning the necessary deviations 
from exponential decay (see e.g. Refs. dHUI) may there- 
fore be confronted with our current theoretical and ex- 
perimental understanding of decaying quantum systems. 
Concerning experimental studies of deviations from the 
conventional exponential decay we notice in particular 
the study of decaying r-leptons the observed de- 
viations in quantum-mechanical tunneling processes [l7| 
and the power-law behavior of the decay at times-scales 
larger than twenty lifetimes in dissolved organic mate- 
rials [l8j | . In nuclear physics the decay of thorium has 
also been s ugg ested as a potential target for large-time 
deviations |19j . 

In the context elementary particle physics gauge in- 
variant definitions of observable quantities are of central 
importance and the definition of a the decay width of 
an unstable particle is highly non-trivial (see e.g. [2(|) 
in this respect. In studies of a possible proton decay 
in Nature the time-dependence of decaying systems may 
also play an important role [2l|. Recently it has also 



been speculated that a non-exponential decay of slowly 
decaying 14 C nuclei may be of practical importance in 
radio- isotope calibration methods [221 ] - 

In the presence of material bodies with non-trivial dis- 
persive properties, it has been observed that the decay 
time of atoms can be changed by many orders of mag- 

which also has been ob- 
|. A natural question to 
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nitude (see e.g. Refs. [2j 
served in the laboratory 
be considered could then be to what extent deviations 
from an exponential decay of atoms or molecules can be 
observed, a subject which has been addressed in great 
detail in the literature (see e.g. Refs. (33l - [37l l). 

In the present paper we will investigate possible de- 
viations from exponential decay as studied in terms of 
an effective two-level system interacting with a continu- 
ous number of degrees of freedom of the electro-magnetic 
field. Analytical as well as numerical methods will be em- 
ployed and will be shown to lead to a consistent picture 
of the small and large time deviations from a constant 
decay rate. 

The paper is organized as follows. In the next section 
we outline the theoretical framework and explicitly verify 
unitarity of the two-level system. A regularized version 
of the integral equation for the decay amplitude is de- 
rived in Section [ill Al and exact short-time and long-time 
expansions are obtained in Section MI Bl Laplace tech- 
niques are employed in Sections 1111 CI and MI Dl in order 
to compare with previous asymptotic expansion results 
as given in the literature and the important role of the 
regularization procedure made used of is emphasized. Fi- 
nal remarks are given in Section IIVI and various techni- 
cal details of some of the calculations are, for the readers 
convenience, summarized in Appendixes IA1 [Bl andlCl 
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II. GENERAL THEORY 

Let us consider a neutral atom at a fixed position ta- 
In order to be specific, we will explicitly consider hyper- 
fine interactions but the results obtained can easily be 
rephrased in terms of clcctric-dipole interactions. The 
magnetic moment of the atom interacts with the quan- 
tized magnetic field via a conventional Zeeman coupling. 
The total, un-renormalized, Hamiltonian has then the 
standard form 

H = E] Tiuj a | a) (a | 

a 

+ J ' d 3 r J dwRwft(r,a;) -f(r,w) + H' , (1) 



where the effective interaction part is 

h' = -EE wi/^-bm 



(2) 



The Hamiltonian H can be regarded as a low-energy ef- 
fective description of a more fundamental and renormal- 
izable theory of electro-magnetic processes, i.e. quantum 
electrodynamics. Here f(r, uS) is an annihilation operator 
for the quantized magnetic field, \a) denotes the atomic 
state and E a is the corresponding energy. We assume 
non-degenerate states, i.e. E a ^ E@ for a / ^. The 
magnetic moment of the atom is fx a p = {<AiAP)i where ft 
is the magnetic moment operator. The magnetic moment 
will typically be of the form /j, a p = gs(a\S/h\(3)eh/2m e , 
where S/h denotes a dimensionless spin operator and 
gs is an appropriate gyro-magnetic factor. The quan- 
tized magnetic field B(r) = B( + '(r) + B^^(r) is ex- 
pressed in terms of B^ + ^(r) = V x A^ + ^(r), where 
B(~)(r) = (B( + )(r))t, and where the vector potential 
is given by 



x G(r,r>').f(r',a/) . (3) 

Here the imaginary part of the complex permittivity is 
e/(r,w) and obeys the Kramer-Kronig dispersion rela- 
tions. The dyadic Green tensor G(r, r',w) is the unique 
solution to the Helmholtz equation 



^ x ^ x G( 



-e(r,u>)G(r,r',u) 



= <5(r-r')l, 



(4) 



where the arrow in ^ denotes a derivation with respect 
to the first argument in the dyadic Greens function. Since 
the Helmholtz equation is a linear differential equation, 
the associated Green's tensor can be written as a sum 
according to 

G(r, r', u) = G°(r, r', u) + G S (r, r', u) , (5) 



where G°(r, r',ui) represents the contribution of the di- 
rect waves from a point-like radiation source in an un- 
bounded medium, which is vacuum in our case, and 
G (r, r',cj) describes the scattering contribution of mul- 
tiple reflection waves from the body under consideration. 
The presence of the vacuum part G°(r,r',w) in Eq.© 
will in general give rise to divergences. A regulariza- 
tion prescription is therefore required. For electric-dipole 
transitions it is well-known that it is sufficient to subtract 
an energy shift, corresponding to the introduction of a 
renormalized mass, which leaves us with a logarithmic de- 
pendence of a cut-off frequency. For magnetic transitions 
this subtraction procedure is, however, not sufficient as 
will we discuss in more detail below. Our strategy is to 
allow for a sufficient number of subtractions to generate 
a logarithmic cut-off dependence also for magnetic tran- 
sitions. 

For reason of simplicity, we now limit our attention 
to a two-level atom approximation, i.e. an atom with an 
excited state and a ground state with the frequency tran- 
sition uja = (E e — E g )/h > 0. We consider the Hamil- 
tonian in Eq.(fT]) and apply the well-known Weisskopf- 
Wigner theory for the transitions e — > g. The solution to 
the time-dependent Schrodinger equation in the rotating- 
wave approximation (RWA), i.e. applying H' « Haw A, 
where 

Hrwa = -\e)(g\fi eg -B^(r A )+h.c. , (6) 
is then (uj g = E g /h and uj e = E e /K) 



m)) = ce(t)e" 



|e) ® |0) 



d 3 r I duj 



3 

E 

m— 1 



x c ff (r,w,m|t)e -i (' u+w «J t |ff)®|r,a;,m) , (7) 

as we have ignored any higher order photon state than 
the 0-photon and 1-photon state. Here |e) and \g) are 
the excited state and ground state for the atom, re- 
spectively. Furthermore, |0) denotes the vacuum of the 
electro-magnetic field and r, oj, m) — /^(r, uj) |0) is a one 
photon state. We will restrict ourselves to the initial con- 
ditions c e (0) = I and c g (r,u;,m\t = 0) = 0. More general 
initial conditions will require a more cumbersome analy- 
sis and will not be discussed in the present paper. The 
probability amplitudes c e (t) and c g (r, ui, m\t) must now 
obey the unitarity condition 



\c (I)' - / d 3 r J duj E Mr »i D\ 



m—l 



(8) 



The probability amplitude for the excited atomic state 
c e (t) is then determined by (see e.g. Refs. 38, 39]) 



dc e (t) 
dt 



dt' K(t-t')c e (t'), 



(9) 
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where the un-regularized kernel K(t) is 



where 



K(t) 



1 

2^ 



due-^-^rtr^u) . (10) 



Here we have made use of the conventional spin-flip decay 
rate for spontaneous emission as given by 

T(r, w) = ^ ^ • Im[ V* x G(r, r, w) x V] • , (11) 

where the arrow in V now denotes a derivation with re- 
spect to the second argument in the dyadic Greens func- 
tion. Below we will also make use of the expression 



c g (r, uj, m\t) = -u) 



i //i e/(r,w) 



71"/) 



^ xG*(r A ,r,Lj))] I(u-u A ,t), 



where we have defined the integral 



J(w,t)= / dt'e lU)t c e (t'). 
Jo 



(12) 



(13) 



The useful identity 
3 f uj 2 

^2 / d3r — £ i( r i UJ ) G 3m(r',r,uj)Gi m (r",r,uj) = 

m—l 

ImG j7 (r',r'», (14) 
and Eq. p^|) now lead to 

/> poo 3 

/ d 3 r / dui \cg{v,uj,m\t)\ 2 = 
1 f°° 

— J dujT(r A ,uj)\I{uj~uj A ,t)\ 2 . (15) 

It now follows from Eqs.© and (fTS"|) that 

r r°° 3 \ 
e {t)\ 2 + d 3 r dujY] \c g (r,w,m\t)\ 2 = 0, 

J J ° m=l / 



dt 



(16) 



i.e. the unitarity condition Eq.© is fulfilled for all times 
as it should. We notice, what may appear to be a triv- 
ial fact, that general time-dependent phase- redefinitions 
of the amplitudes will not change the unitary condition 
Eq.©. This circumstance will, nevertheless, be useful 
below in order to circumvent divergent vacuum fluctua- 
tions of the theory. 

Eq.©, a well-known Volterra integral equation of sec- 
ond kind, may be integrated with respect to time. The 
result is then 



c e (t) = l+ dt' K(t - t') c e (t') , (17) 
Jo 



n(t) = / dt'K(t'), 



is the time-integrated kernel, i.e. 



«(*) 



1 

2^ 



du ■ 



i(u - uj a ) 



(18) 



T(t a ,uj). (19) 



In passing, we mention that a rate similar to Eq. (jlip 
may be derived for electric-dipole transitions (see e.g. 
Refs. [3ll - l4ll |). In that case, the spontaneous decay rate 
r(r, w) should be replaced by 



d eff -Im[G(r,r,u;)] -d ge , (20) 



where d eg is the electric dipole moment for the tran- 
sition e — > g. By comparing Eqs.flTTj) and (|2T))) it 
is now straightforward to translate between magnctic- 
dipole and electric-dipole transitions which will be made 
use of below. When referring to electric-dipole transi- 
tions we therefore assume that one make use of the ap- 
propriate decay rate, i.e. Eq. (l2"0)l . 



III. VACUUM 



Exact 



The dyadic Green tensor for vacuum, G (r,r,q/), is 
^iven by a well-known expression (see e.g. Refs438, 39]) 



Im[^ x G°(r,r,w) x V] 



6ir c A 



1 
1 
1 



. (21) 



independent of the position r as it should due to the 
translational invariance of the vacuum quantum state. 
In this case the kernel in Eq. (TT9l reduces to the vacuum 
kernel K n (t) and reads 



Tn f Uc UJ e -i(t»-WA)t _ l 



2tt 



wa i(u)-u) A ) 



(22) 



with n = 3, and where we have introduced a cut-off fre- 
quency ui c in order to make the integral finite. The decay 
rate of magnetic spin-flip transition for a two-level atom 
with no angular momentum and negligible nuclear mo- 
ment in free space is then (see e.g. Ref . [24| ) 



with 



fJ-o 



37T Ti 



(23) 



(24) 
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and kA = ^a/c is the wave number in vacuum. Here 
we have introduced the dimensionless spin factor S 2 = 
S 2 + S 2 + S 2 , where Sj = (g\Sj/h\e) is the dimensionless 
matrix element component of the electron spin operators 
Sj corresponding to the transition |e) — > \g), with j — 
x,y,z. Furthermore, /is = eh/2m e is the conventional 
Bohr magneton, gs ~ 2 is the electron spin gs factor. 

Clearly, the time-dependent kernel in Eq. (j2"2")l is diver- 
gent as uj c — > oo and it is not entirely clear to us how 
to make sense of this kernel for all times t. Since, as we 
have mentioned above, magnetic-dipole transitions are 
analogues to electric-dipole transitions, we therefore ap- 
peal to Bethe's mass-renormalization [42j procedure as 
far as dealing with divergences are concerned. At large t 
we recall that one then e.g. can make use of the distri- 
butional identity 



-i(u-u A )t _ i / i 
lim = P 

t-¥oo ui — UJA \ UA — W 



ittS(uj — uja) ■ 



(25) 

where one from the principal part P(l/(u> A — oS)) then 
adds a term P(1/uj) which, when summed over all possi- 
ble final states, corresponds to Bethe's introduction of a 
rcnormalized mass. In the electric-dipole case, the finite- 
time kernel K n (t) will then be regularized in such a man- 
ner that, for large t, energy shifts are reduced to at most 
a logarithmic dependence of the cut-off frequency oj c . As 
we will see below in Section IIIIDI this regularization can 
be carried out in terms of frequency shift and a con- 
ventional mass renormalization for all times with a mass 
counter term. 

In the present case of magnetic transitions, which cor- 
responds to an interaction which is not directly renor- 
malizable in the same manner as for electric-dipole tran- 
sitions (see e.g. Ref. (43|) we, nevertheless, proceed in a 
manner which treat these different transitions in an equal 
manner as far as divergences are concerned. We therefore 
subtract the second order expansion of the denominator 
in Eq. (|19p for u) 3> uja in such a way that at most a loga- 
rithmic dependence of the cut-off frequency ui c remains, 
i.e. 



f 



1 



■ UJA 



uj — UJA 



1 



UJA_ 

,.,2 



(26) 



which leaves us with the following regularized version of 
the kernel Eq. 



4(t) = g 

Z7T 



duj 



1 



i(uj ~ uja) 



(27) 



The kernel Kq (t) is therefore obtained from Ko(t) in 
Eq. ([22"T) using the rule r(r,w) — > T(t,uj A ). With regard 
to the amplitude c g (r,ui,m\t) we observe that unitarity 
prevails if we in Eq. (fT5"|) also make use of the same pre- 
scription. The fact that unitarity is preserved is our pri- 
mary motivation for the introduction of the regulariza- 
tion procedure above. Below we will, however, investi- 
gate the effect of keeping a dependence uj /uja instead of 



(uj/uja) 3 in the kernel Eq. (l22l) . This would be in line with 
the standard electric-dipole transition considerations as 
used in e.g. Refs. |33l - l37j |. As we will verify below, the 
choice of regularization procedure will effect the large 
time behavior of the decaying system. 

In Eq. ([27[) we make use of a cut-off frequency ui c in 
order to make the frequency integral finite. Since our 
calculation is non-relativistic we can, e.g., identify this 
cut-off with with m e c 2 /Ti, but one may also regard this as 
a free parameter when one e.g. considers artificial atoms 
(Ref. [4[ and references cited therein). We therefore in- 
troduce a dimensionless parameter cut-off parameter A 
as defined by 



A 



OJc 
UJ A ' 



(28) 



The subtraction procedure of Ko(t) above enables us to 
extract the leading logarithmic dependence of the cut-off 
parameter A , and, in the end, we have therefore replaced 
the vacuum kernel Ko(t) with the following regularized 
kernel K$(i) 



27r J-UJA X 

To [ {A - 1)UA , cos{xt) , To 



2tt 



dx ■ 



i^ln(A-l).(29) 



The imaginary logarithmic term in this equation corre- 
sponds to an induced Lamb shift due to vacuum fluctua- 
tions since it can be removed by introducing an angular 
frequency shift, i.e. 



ua = u A - — hi(A - 1) 
and applying the transformation 



c e {t) c e (t) 



c e {t) exp ( - ln(A 

Z7T 



(30) 



1) , (31) 



as can be seen by making use of Eq. (|9]) . We also find it 
convenient to define the kernel 



«?(*) = - IT [ Si P - 1)wa*) + Si[Q A t] 

-t^ (ci[(A-l)a}^]-Ci[o} A t]^ , (32) 

where Ci(x) and Si(x) are the standard Cosine and 
Sine integral, respectively (see e.g. Ref. [441] ) and A = 
Ku)a/&A — A in view of the fact that A is considered to 
be a large cut-off parameter. The probability amplitude 
c e (t) is now a solution to the integral equation Eq. (|17[) 
provided the kernel K${t) is replaced by h^it), i.e. we 
consider 



o 



c e (t) = l+ dt' Ro(t-t')c e (t') 



(33) 
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The solutions of the corresponding integral equation can- 
not be obtained in closed form and therefore we will re- 
sort to a numerical treatment. The numerical results as 
presented in the present paper involves algorithms with 
arbitrary numerical precision. 



B. Approximations 

It is convenient to define a dimensionless time param- 
eter r = Tot as well as the parameters bA = loa/Tq and 
b A = ma/To, i-e. 

~ bA = bA -^^. (34) 

Z7T 

Let us now consider small time scales such that (A — 
l)bAT <C 1 and b A r <C 1. Such small times are basically 
of academic interest as t <C h/m e c 2 — 1.29-10 _21 s unless, 
as suggested above, a different scale can be provided by 
the use of e.g. artificial atoms. Virtual particles may 
then be created and our theory is not strictly valid. The 
leading series expansion of the kernel as given by 

Eq.(|32[). for such small time parameters t is: 



0.4 



0.2 



I , , 1 , , 1 , , , 1 

0.4 0.8 1.2 1.6 2 

T 

FIG. 1: The probability |c e (V)| 2 = |c c (r)| 2 as a function of 
t = Tot. The solid curve corresponds to the exact numeri- 
cal solution, i.e. solution of Eq. (|33[) with the kernel k^(t) as 
given by Eq. (l32[) in the main text for a two-level system in vac- 
uum. The upper (dash-dotted) curve corresponds to the small 
time, \>at <C 1, expansion Eq. (|39p . The lower (dotted) curve 
corresponds to the exponential decay exp(— r) with expected 
deviations at small r. The values of the relevant parameters 
are b A = u)a/T = 10 and A = 1000. 



~4(r) « - ^ A6 A r - i ^ ln(A - 1) , (35) 

using A w A. Apart from a sign, the imaginary part 
of this kernel is the same as the last imaginary part of 
Eq. (P^)l . It is therefore convenient to invert the transfor- 
mation Eq. (|3ip and consider the probability amplitude 
c e (r). Substituting the kernel Kq(t) into Eq. (fT7|) , we 
obtain 

c e (r)^l-^-(b A rf, (36) 

using the approximation c e (t') ~ 1 in Eq. (|17[) since time 
parameter r is small. We observe that in Eq. (|3"6")) the 
parameter bA enters and not bA- In passing, we also 
mention that an expansion for small times is also given 
in Ref-HH (their Eq.(3.21)) in the case of electric-dipole 
transitions, with the result 

c e {j)^l-^-{b A rf. (37) 

The difference between equations (|36j) and (|37j) . in which 
case, as mentioned above, the factor (lu/uja) 3 in Eq. (p?2")) 
is replaced by uj /loa , illustrates the relevance of the sub- 
traction procedure. Below we will see that this effect 
is even more pronounced at large vales of r. The dis- 
crepancy between these last two expressions is easily ex- 
plained by the fact that different regularization proce- 
dures have been applied. 

Let us still consider small times <C 1 but such 
that (A — l)b A T 3> 1. Such time scales are of relevance 



l 



0.98 



0.96 



0.94 




0.92 



0.901 ' ' ' ' ' ' ' ' ' 1 

0.02 0.04 0.06 0.08 0.1 

T 

FIG. 2: The probability |c e (i~)| 2 = |c e (r)| 2 as a function of 
r = Tot as in Fig. [T] but for a smaller r interval. 



when one e.g. study the Zeno effect for quantum systems 
(see e.g. Refs.fl5l. l45j|). For such intermediate times, the 
kernel in Eq.(|32|) is reduced to 

~4(r) « - ^ + i g ( 7£ + \n(b A r) ) , (38) 

where je ~ 0.577 216 is Euler's constant. Substituting 
Eq. (l55|) into the Eq. dB^l) . expressed in terms of c e (r) and 
using c e (r') w 1, we may then carry out the time inte- 
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FIG. 3: The probability |c e (r)| 2 = |c e (-r)| 2 as a function of 
t = Tot. The solid curve corresponds to the exact numerical 
solution, i.e. solution of Eg. (|33fl with the kernel Ko(t) as 
in Eg. (132fl for a two-level system in vacuum. The almost 
overlapping dotted curve corresponds to Eq. (|41[) . The values 
of the relevant parameters are b A = loa/^o = 10 and A = 
1000. The inset shows a limited part of the larger figure and 
illustrates the accuracy of the asymptotic form Eq.(|41l). 



gration, with the result 



Ce(r) 



r 
4 



2tt 



j e + ln(b A T) - 1 



(39) 



The corresponding probability |c e (r)| 2 = |c e (r)| 2 is illus- 
trated in FigsfJ] and [5] (dash-dotted line) together with 
the exact numerical solution, i.e. the numerical solution 
of Eq.® with the exact kernel Ea.([3"2| (solid line). 

Finally, for large times as compared to the shifted 
atomic frequency transition, i.e. b A r 1, and also 
(A - 1)Iat > 1, the kernel Ea.([3"2"]l is, to first leading 
order 



U V I r, 1 1_ 7 



, -L 



r n e ibA7 



2 2tt b A r 



(40) 



One may then show that (see next section and Appendix 



--t/2 _ 



2niT[ b A - \n{b A T)/2-K 



, (41) 



provided b A is sufficiently large, and where we have made 
use of A « A. This solution includes the well-known ex- 
ponential decay as well as a correction term which dom- 
inates |c e (T) | 2 for large values of r. Eq. (|4"Tj) is plotted 
in Fig. [3] (dotted line) together with the exact numerical 
solution (solid line). The last term in Eq. (|4Tj) will dom- 
inate over the exponential for times r > t* , where r* is 



determined by the transcendental equation 

1 



,-r*/2 



2-kt* b A - \n(b A T*)/2ir 



(42) 



valid for a sufficiently large value of b A . For b A = 10 
and A = 1000 as in Fig. [31 the solution of this equa- 
tion is, e.g., r* « 18.5. The t\\\'(J) A t) dependence in 
Eqs. (j4TT) and (|42"T) becomes more important than the r 
dependence only for very large time-scales r > Ti n , where 
Ti„ = exp(27r6 y !i)/6 y !i, in which case |c e (r)| 2 becomes ex- 
ponentially close to zero with increasing value of 6,4- 

C. Asymptotic Expansion 



Laplace transform techniques have been used to inves- 
tigate possible deviations from exponential decay as in 
Refs. [33L l35l. l36l] . As shown in particular by Seke and 
Herfort [35|, [3fJ| , careful considerations of the analytical 
properties of the decay amplitude are required in order 
to extract the asymptotic behavior of c e (t). The result 
of Refs.JH, |H, Hg| may seem to be partly contradictory. 
It is therefore of some interest to investigate in detail 
in what manner the asymptotic behavior according to 
Eq. (jTl|) is obtained by making use of Laplace transform 
techniques. As we will see, it then becomes apparent that 
the various asymptotic forms of the decay amplitude are 
all valid but depends on the time-scale and other physical 
parameters at hand. 

In order to obtain the corresponding asymptotic ex- 
pansion, as r becomes large, we therefore find it conve- 
nient to consider the Laplace transform of the differen- 
tial equation for the rescaled amplitude c e (t) as defined 
in Eq.dnj), i.e. 



dc e (t) 
dt 



dt' Kg(t-t')d e (t') 



igln(A-l)ce(i), 



(43) 



where the regularized kernel K^{t) is 



2tt 



uj a 



(44) 



with, for the moment, n = 0. The Laplace transform 
c e (s) is then given by 



c e (s) 



1 



s - Kg(s) + iTo ln(A - 1)/2tt 
where we have identified 

- v r dL ° — ^— = 

Air J q u! — lo a — is 
^fln(l--^' 

27T V IS + uja 



(45) 



(46) 
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in terms of the principal branch of the natural loga- 
rithm ln-function. The inverse Laplace transform in 
terms of the Bromwich integral and the dimensionless 
time parameter r can then formally, after a suitable Wick 
(u = b A + is/Tq) rotation, be written in the form 



gib at 

2m 



due' 



00+27 



u-b A + ^-ln(l-^) 
Zir u 



where we have defined 



r 



(47) 



(48) 



and where the positive real number 7 is chosen in such 
a way that the possible singularities of the integrand are 
above the integration contour. 

We now follow the methods of, in particular, Refs.(35l. 
|36| . where more details can be found, and use Eq. (|4"7|) 
in order to extract the large r, i.e. u — > 0, asymptotic 
expansion for c e (r). As in Refs.[35l we identify the 
functions 



M (u) 



A c )-log(u)], (49) 



and 



Mt{u) = M (u)+i, 



(50) 



where the log-function now stands for the multi-valued 
natural logarithm function. The appropriate Riemann 
surface for the function log(u — A c ) — log(u) has been 
obtained in Ref . (35j| . Here one joins a second Riemann 
sheet at the branch-cut along the real M-axis from u = 
to u = A c . 

The integration contour in Eq. (|4"7|) can then be de- 
formed into two curves of integration that runs paral- 
lel to the Im(u)-axis. One of these curves runs from 
— ioo+A c -|-e , above and around the branch point u = A c 
on the first Riemann sheet, and then towards — ioo+A c — e 
on the second Riemann sheet, where e is a small and posi- 
tive real number. Similarly, the other curve of integration 
runs from — ioo + e on the second Riemann sheet, above 
and around the branch point u = on the first Riemann 
sheet, and then towards — ioo — e on the first Riemann 
sheet. 

The function M\ (u) is obtained from Mo (u) when pass- 
ing to the second sheet of the Riemann surface con- 
struction. It can be shown that each of the functions 
Mq(u) and Mi(u) have no poles and onl y o ne zero on 
each of the Riemann sheets considered [35l . [36| . On 
the first Riemann sheet, where the original functions in 
Eq. (|4"T)) are defined, there is, in addition, a zero of Mo (it) 
along the real axis at uq « A c (l + exp[— 27rA c ]). The 
corresponding residue is exponentially suppressed, i.e. 
Z = l/(dM (u)/du) u=Uo « 27rA c exp(-27rA c ), apart 
from a phase factor, and the corresponding pole contri- 
bution has therefore been neglected. 



In the deformation of the integration contour in 
Eq. (|4"T|) , and on the second Riemann sheet, one now 
encounters the conventional Wigner-Weisskopf like pole- 
contribution at u = u\ « b A — i/2, provided b A 3> 1 due 
to Eq.([50]). Since dMi(u)/du = 1 + A c /2ttu(u - A c ) « 
1 — l/b A w 1 for u = u\ provided b A 3> 1, this pole 
leads to a residue Z\ = l/(dMi(u)/du) u — Ul « 1. This 
Wigner-Weisskopf like pole-contribution therefore con- 
tributes to the amplitude c e (r) with exp(— t/2), i.e. the 
expected and conventional amplitude describing expo- 
nential decay. Including the two integrals emerging from 
the branch points mentioned above we then, finally, ob- 
tain that the result 



c e (r) 



-T/2 



a^AT 



-h 



3 i6 A (l-A)r 



27TT 2tTT 

where we have defined the integrals 



h = 



dse 



1 



1 



and 



Mo(-is/r) Mi{-is/r) 



h = 



(51) 



(52) 



dse 



Mi (A c - is/r) Mo (A c - is/r) 



(53) 



similar to the results of Refs.fl^, [36|. In Appendix [B] we 
now show that, as t > 1, 



h 



(b A -Hb A T)/2^ 



(54) 



provided \b A - ln(A c r)/27r| w \b A - \n(b A T)/2n\ > 1. In 
the same manner it also follows that 



(55) 



(Ac - b A - ln(A c r)/27r) 



provided |A C — b A — ln(A c T)/27r| 3> 1. As long as A c b A 
it is now clear that I\ will dominate over li- Due to 
Eq. ([54|) we now recover the result as given in Eq.pTj). 



It now also follows that l\ w l/b A unless r is exponen- 
tially large, i.e. r 3> e 2vbA /b A , in which case the decay 
amplitude c e (r) is exponentially small. 

When compared to the analysis of Refs.(35l. |3(| we re- 
mark that our approach to the asymptotic behavior of in- 
tegrals like Eqs. (|52|) and (IB1[) has a much broader range 
of applicability. In particular we notice that, as long as 
t < e 2rrbA /b A ), the h integral leads to a 0(1/ t) behav- 
ior of the decay amplitude c e (r). It is only in theoppo- 
site limit that the proposed C(l/rln 2 T) of Refs.[35l l36j 
emerges, i.e. for r 3> e 2lxhA jb A . In this case the decay 
amplitude, however, becomes exponentially small for b A 
sufficiently large. We also find it important to express 
the asymptotic behavior in terms of dimensionless quan- 
tities like r and not in terms of the time variable t with 
a physical dimension as in Refs. [35l [36| . 
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FIG. 4: The probability |c e (r)| 2 = |c e (r)| 2 as a function of 
r = Tot. The upper curve corresponds to the exact numerical 
solution, i.e. solution of Eq. ()33p with the kernel ita(t) as in 
Eq. (|32[) for a two- level system in vacuum which, within the 
numerical accuracy of the figure, overlaps with the asymp- 
totic formula Eg. 1)411) . The lower curve corresponds to the 
main result in Ref . [331] . i.e. Eg. (l73|) . The asymptotic form ac- 
cording Eq. (|75|) is in excellent agreement with the numerical 
evaluation of Eq. ()73p . The values of the relevant parameters 
are 6 a = 10 and A = 1000. As a reference, the straight line 
in the Figure corresponds to a pure exponential decay, i.e. 
|c e (r)| 2 = e-. 



1M 



D. Comparison with Refs.[3 



As alluded to above, in many cases of the study of 
deviations from the exponential decay of atomic spon- 
taneous emission processes, various approximations have 
been applied in order to obtain explicit and analytical ex- 
pressions as, e.g., in Refs.f33l.l35l.l36j and references cited 
therein. In these references one consider electric- dipole 
transitions and, in the integral equation Eq. (ll7[) for the 
vacuum case, the kernel k{t) should therefore be replaced 
by kernel K\{t) = «i(t) as defined by Eq. (f2"2")) . i.e. with a 
factor lu/lja in the integrand rather than (lo/lua) 3 - Con- 
trary to the regularization procedure employed in Section 
IIII1 this factor oj/loa is, however, kept untouched and a 
suitable subtraction is carried out at a later stage in the 
analysis. As we now will verify the actual subtraction 
procedure in order to make the corresponding kernel fi- 
nite will effect the large time deviation from exponential 
decay. If the cut-off frequency uj c again is regarded to be 
a finite physical quantity, we notice that the regularized 
kernel K\{f) may be written in terms of k^(t) as given 
by Eq.([27D: 



2lTbAT 



«f(r) 

-i(A-l)b A T _ 



„ibAT 



i^A. (56) 



linear in the dimensionless cut-off frequency rather than 
the logarithmic dependence of Kq(t) as in Eq. (j2l?)) . We 
observe that the linear A-term may now be absorbed in 
an energy shift such that the kernel Ki (r) is replaced 
according to: 

«? (r) -> «f (r) = 

K R( T ) + _iiL_ ( e -i(A-i)b' A r _ e «i 
2tt W a t V 

where b' A = uj' a /Tq, and 



(57) 



U) A = UA 



£0 

2tt' 



A 



(58) 



and where also the w^-dependence in k,q(t) is replaced 
by w'j . Such a frequency shift preserves unitarity. As in 
Refs.jH, H^, |36j], the linear A-dependence can therefore 
be removed altogether by applying Bethe's mass renor- 
malization [42j , i.e. by including a suitable mass counter- 
term in the Hamiltonian. After this removal of the linear 
A-dcpcndence and by, in addition, applying the transfor- 
mation in Eg. pip, the kernel k^(t) is transformed into 
the kernel £f (t) according to: 

«f (r) «f (r) = 

jj*( r ) + ( e -*(A-i)5Ar _ e fb A r\ (5Q) 

The asymptotic form of the kernel K± (r) is now such that 
the sub-leading contribution from the asymptotic from of 
the kernel Rq(t), as given by Eq. (|4T)l) . is cancelled and 
therefore 



2 27T & A r 



(60) 



for sufficiently large r. Hence, the methods of Appendix 
lAl would then lead to 



C e (T) 



-r/2 



,-i(A-l)b A T 

2iriTkb 2 A 



(61) 



apart from possible 1/r 2 corrections provided &a is suf- 
ficiently large. 

As in Section [ill CI a more precise asymptotic expan- 
sion as r becomes large, and which also takes the Lamb 
shift Eq. (j3"0")) more carefully into account, can be obtained 
by a consideration of a Bromwich integral. The Laplace 
transform of the differential equation for the rescaled am- 
plitude c e (t) reads 



dt 



dt' K?{t-t')c e {t') 





r 



$TYl C 

1— ln(A — l)ce(t) — i — — — c e (t) , (62) 
where the regularized kernel K^(t) is 



Km = 



2ir 



du> — e 



Due to this identity we realize that the kernel t^(r) is 



(63) 
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and where the unitarity preserving transformation as de- 
fined by Eq. (j3"Tj) has been applied. In Eq. (j6"2")l we have 
also included the contribution due to the additional mass 
counter-term 5m c 2 /% in the Hamiltonian which has the 
role of cancelling the linear A-dependence in the kernel 
Eq. (|63| . After a tedious but straightforward analysis, 
following the same Laplace transformation techniques as 
in Refs. [35l |36|| . the large r expansion of c e (r) can now be 
obtained. Similar to the definitions of Mq(u) and Mi (it) 
in Section [III CI we therefore introduce the functions 

N {u) =u-b A + ( log(u - Ac) - log u) , (64) 

and 

Nx{u) = N {u) + i— , (65) 

OA 

where the linear A-dependence in Eq.(|62j) now has ex- 
plicitly been removed by the mass renormalization pro- 
cedure. The solution to the Eq. (R>2")) can then be written 
in the form 



c e (r) = 



ib A T r — ca+il 



2ni 



due 



1 



oc+27 



N (u) 



(66) 



By deforming the integration contour in Eg. (1661) in the 
same manner as in the derivation of Eq. (|51[) . one then 
obtains the result 

p ib A T ib A (l-K)r 

~^) = e-^-—J l --^—J 2 , (67) 
where we have defined the integrals 



oo / I 

— s ( 



Jl Jo dSe ~ S {N (^s/r) mi-is/r) 



. (68) 



and 



J 2 = 



dse 



N^Ac-is/r) N (Ac-is/T) 



(69) 



We have therefore actually reproduced the results of 
Refs. [35, 36], but now expressed in terms of dimensionless 
parameters. 

As expanded upon in the Appendix O the integrals 
Ji and J2 can now be analyzed in manner similar to the 
integrals I\ and I\ of Section ITlI CI . If r 3> 1 and A c 3> 1 
and, in addition, fe^r > 1 + ln(A c )/27r&yi , we then obtain 



For the integral J 2 it also follows that, for t>1, 

1 



■h 



iA(bA — ln(A c r)/27r) 2 



(70) 



(71) 



provided \oa — ln(A c T)/27r| ^> 1. The insertion of these 
asymptotic expansions into Eq. (|67j) therefore leads to 



c e (r) 

e -i(A-\)b A T 1 



-t/2 



J.b A T 



1 1 



2nb 3 A r 2 



2niA t (b A ~ln(A c T)/27r) 2 



(72) 



The last term in Eg. ([72]) with the (b A ~ ln(A c t)/2tt)- 
dependence has been obtained by using the technique 
of the Appendix [B] in order to see how the claimed 
rln 2 t dependence of Ref.^S, ,36] actually emerges. Here 
we again stress that our results are expressed in terms 
of dimensionless quantities in contrast to the results of 

Refs. m, m. 

At this point we now quote the following explicit, but 
approximative, expression for the probability amplitude 
c 6 (t) result of Knight and Milonni in Ref.[33j: 



Ce(r) = e - T / 2+l&bT 

O 

A X bA t 



= C: 

I dx e 



e ib A T 
2lTi 



1 



X A + 1 2b A 



2tt b A 



[In a; - 



where 8b = Suj/Tq, and 



Su = ^ ln(A - 1) , 



(73) 



(74) 



is the angular frequency shift of the excited state \e). 
This expression, which we have verified [46| . is obtained 
by deforming integration contours of the Bromwich inte- 
gral by only considering the analytical properties of the 
natural logarithm ln-function and not taking into account 
the additional analytical conditions imposed by the ap- 
pearance of the functions N Q (u) and N\ (u) as discussed 
above. For lar ge t imes, i.e. &at ~3> 1, this expression is 
reduced to 111, Hal 



Ce(r) 



D -r/2+i8bT 



1 



1 



Jb A T 



2-Kb A (b A T) 2 



(75) 



We notice that up to a phase redefinition, this equa- 
tion and Eq. ([72|) only differ by that last Seke-Hcrfort 
term of Eq. (PT2")) . It is only in the limit A — > 00 that 
Eq. (|72|) leads to 1/r 2 asymptotic behavior of Knight and 
Milonni. It is now also clear that it is only for time-scales 
t > e 2rrbA /A c that the 1/r In 2 r asymptotic behavior of 
Seke-Herfort emerges. Such time-scales are typically ex- 
ponentially large leading to an exceedingly small decay 
amplitude c e (r) unless A c is very large. But as we have 
seen above, in such a case the 1/r 2 asymptotic behavior 
of Knight and Milonni begins to dominate. We also re- 
alize that Eq. (f72")) leads to a dominant 1/r behavior for 
r < e 2lTbA /A c unless A c again is not arbitrarily large. 
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FIG. 5: The probability |c e (r)| 2 = |c e (r)| 2 as a function of 
r = Tot. The upper curve corresponds to the exact numerical 
solution of the Volterra equation as considered in Ref . [35l. |36| . 
i.e. of Eg. ((33)) with the kernel hf(r) of Eq. ([59)) for a two- 
level system in vacuum. For large values of r the solution 
is dominated by the 1/r asymptotic behavior according to 
Eq. ([72]) . The lower curve corresponds to Eqs. (f73|) and ([75]) . 
which overlaps within the numerical accuracy of the figure, 
with a characteristic 1/r 2 asymptotic behavior for c e (r). The 
values of the relevant parameters are b A = 1000 and A = 1000. 



By a comparison of Eqs. ([41]) and (|72[) we also realize, 
due to the presence of the factor 1/ A in the Seke-Hcrfort 
term in Eq. (]72p , that the regularization procedure actu- 
ally effects the large-time behavior of the decay ampli- 
tude c e (r). 

We illustrate some of these features in Fig. @] and 
Fig. [5] For comparison, Eqs. ([73|) and ([75]) are plotted 
in Fig. 0] together with the exact numerical solution. In 
Fig. [5] we have plotted the solution of Eq. ([3"3"]) with the 
kernel hf (r) of Eq.® together with Eqs. ([73]) and ([75]). 
As long as the Seke-Herfort term in Eq. ([72"]) dominates for 
sufficiently large values of r, this contribution will dom- 
inate the conventional exponential decay term if r > r„, 
where r* is determined by the transcendental equation 



-r„/2 



1 



27tAt* ( b A - ln(A c r*)/27r 



(76) 



With the choice of parameters as in Fig. [5] and r > r* ~ 
53, one finds that Eq.([72]) in this case leads to 1/r scaling 
rather than the 1/r 2 asymptotic scaling for the decay 
amplitude c e (r) according to Eq.([73]). 



IV. FINAL REMARKS 



We have seen that the procedure of regularization 
the decay amplitude c e (r) actually influences the devi- 



ation from exponential decay at large values of the time- 
parameter r. It would, of course, be of interest if one e.g. 
in terms of artificial atoms could experimentally inves- 
tigate deviations from exponential decay at large times 
and thereby infer a possible presence of the Seke-Herfort 
term of Eq. ([72]) or the related 1 /r behavior obtained in 
Eq. (l41[) . As alluded to in the introductory remarks and 
in the case of atom chips, the presence of material bod- 
ies can drastically change the characteristic lifetime (see 
e.g. Refs. [23l - [27l ] and one could then perhaps expect to 
encounter deviations from exponential decay. With the 
parameters of Refs-dKH,^, i.e. S 2 = 1/8, g s « 2 and 
uj a /2it = 560 kHz, we have that b A « 3 • 10 32 . Using the 
cut-off A according to Eq.([28l), i-e. A ss 2 • 10 14 one finds 
that b A ~ b A and that r* w 315 and r* w 380 accord- 
ing to Eqs. (HU and ([76"]) respectively. Since our analysis 
involves the dimensionless parameter r this means that 
even though the effective decay-rate T can be changed 
by , say, twenty orders of magnitude (see e.g. Refs. (23l — 
|27|), one still would have to consider several hundred 
decay times in order to see deviations from exponential 
decay. This means that for atom chips the conventional 
exponential decay description applies with an exceedingly 
high accuracy. As we, however, have remarked in the in- 
troduction, experimental studies have e.g. shown that 
artificial atoms can lead to a Lamb shift of the order of a 
few per cent of the typical emission line [J] . We now esti- 
mate the corresponding relevant parameters according to 
Ref.[3|. The natural frequency is u> a /2tt rj O(GHz) and 
the effective rate Tq is estimated from the Lorentzian fit 
to the experimental data to be Tq 0.5 • 1{)~ 2 uj a /2 , k. 
With a Lamb-shift of the order of one per cent of u A , 
we then have that To ln(A — l)/27r ks lO -2 ^, i.e. the 
effective cut-off A is then given by A w e &v . One then 
finds somewhat more optimistic numbers r* « 40 and 
r» rs 200 according to Eqs. (fl*2"]) and ([76]) . respectively, 
which may open the door for observing deviations from 
exponential decay in system composed of artificial atoms. 
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Appendix A: Large time expansion 

In section [ill B[ the kernel Kq(t) was considered for 
three different regimes: small, intermediate and large 
times r for a given b A and given A. The integral in 
the Volterra equation Eq. (l33|) may then trivially be split 
into these three intervals, where one interval relieves the 
other. For large times t> 1 and, in addition b A ^> 1, 
two of these integrals are negligible. The Volterra equa- 
tion in Eq. (|33[) is then left with an integral corresponding 
to the kernel in Eq. (|40|) . This integral may be computed 
analytically: 



r-l/b A 

c e (r) — 1 rj / dx 



1 1 e^Afr-a;) 

2 + 2^b A (T-x) 



-x/2 



= -1 + e -/ 2 +i/2bA _ _J_ El[ _( 1/2 + ib A )T]e~ T ' 2 
2irbA 

+ ^- El [-(l/2b A + l )}e-^\ 

(Al) 

where E\{z) is the exponential integral for a complex 
argument z. As r > 1 and b A 3> 1, the latter expres- 
sion may now be simplified by making use of the leading 
asymptotic expansion of the E\ functions in Eq. (|A1|) for 
large values of the parameter r: 



c e (r) 



,-t/2 



o ib AT 



1 



2ni b\rl- i/2b A 



(A2) 



i.e. Eq. (|4Tj) if b A 3> 1. As alluded to in the main 
text Eq. (|A2j) remains valid for large values of r provided 
6,4 hi(b A T)/2ir, a condition that emerges from further 
iterations of the Volterra equation in Eq. ([3"3"l) . As dis- 
cussed in Appendix [Bj a more precise asymptotic form 
than the one given by Eq. (|A2p can be obtained by mak- 
ing use of Laplace transform techniques. 



Appendix B: Asymptotic Expansion of Ii and I2 

Here we first consider the integral I\ as defined in 
Eq. (f52")l in Section UlI CI which we rewrite in the following 
convenient form 



h = I dse 
lo 



1 



-is/r - b A + ln(l + A c r/is)/27r 
1 



-is/r -b A + ln(l + A c T/is)/2n + i 
i.e., as r becomes large, 

1 



(Bl) 



L w - 



a + lns/27r + ^'/4 
1 



a + lns/27r - 3i/4 



(B2) 



where we have defined 

ln(A c r) 



b A 



2tt 



b A 



2ir 



(B3) 



if A « A. One now finds, by a combination of analytical 
and numerical methods, that the real part of I\ can be 
neglected, provided a is big enough, i.e. \a\ > 10, and 
hence l\ « if a 2 . By making use of the same reasoning 
as above one can now verify that 1 2 ~ i/a 2 , where now 
a = A c — b A — ln(A c r)/27r with A c = lj c /To. Therefore 
the contribution due to I2 can be neglected in comparison 
with Ji, for any reasonable range of r, provided that 
A c > b A . 



Appendix C: Asymptotic Expansion of Ji and J2 

Here we first consider the integral J\ as defined in 
Eq. in Section IIIIDI for large values of r where we, 
due to the presence of the factor e~ s in the integrand of 
J\i make use of the approximations 



No(-is/r) w -b A \1 + i—S ) , 
with (5=1+ ln(A c )/27r6A, as well as 



N^-is/r) w -b A ( 1 + ) + -1- 

b A r J b A r 



(CI) 



(C2) 



as r becomes large. A series expansion in s/b A r then 
leads to the convergent integral 



°4 T Jo 



dse 



(1 + t^/6 A r) 2 



(C3) 



The effective expansion time-parameter is then actually 
b A r instead of r. We therefore find that J\ w l/b A r 
for sufficiently large values of the parameter b A r , i.e. 
b A r > 1 + ln(A c )/27r6A- In the opposite limit J\ will be 
suppressed by factors like \n(5 /b A r) / (5 /b A r) 2 in addition 
to the l/b A r behavior. More care has to be taken when 
considering the integral J2. For large values of r, we then 
make use of the expansions 



Nq{A c - is/r) « A 



In s 
27 



with a as in Eq. (|B3l) . as well as 



3 In s 

N 1 (-is/T)^A[a+-i + — 



(C4) 



(C5) 



Here we have made use of the simple fact that A c = Ab A . 
As long as \a\ is big enough it can now be argued that 
we can neglect lns/27r term in the N and Ni functions 
above. This is so since for large s > e 27r ' a ' the integrand 
of J2 is exponentially suppressed while for sufficiently 
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small s < e -27r|a| ^here will only be an exponentially 
small range of integration. Proceeding in this manner 
and by series expansions in —i/Aa and 3i/4a, we then 
find that the dominant contribution to J 2 is imaginary 



and that J 2 ~ 1/iAa 2 . By numerical methods we find 
that this is indeed an excellent approximation if \a\ > 10 
when evaluating \c e (t)\ 2 . 
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